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Abstract: SU(2) gauge theory with two fermions transforming under the adjoint rep- 
resentation may appear conformal or almost conformal in the infrared, and is one of the 
candidate theories for building models for technicolor. Early lattice Monte Carlo studies 
of this model have used unimproved Wilson fermion formulation, which can be expected 
to have large lattice cutoff effects. In this paper we present the calculation of the 0{a) 
improved lattice Wilson-clover action of the theory. The Sheikholeslami-Wohlert coefficient 
has been determined non-perturbatively, and various boundary improvement terms, needed 
for the Schrodinger functional formalism, have been calculated in perturbation theory. For 
comparison, we have also determined the improvement coefficients for SU(2) gauge the- 
ory with two fundamental representation fermions. The calculation paves way for more 
accurate lattice Monte Carlo analyses of the theory in the future. 



Keywords: 



Lattice field theory, Conformal field theory. 



Contents 



[T]. Introduction |l| 
|2]. Lattice formulation: the model and O(a) improvement ||] 

||. Perturbat ive analysis of the boundary improvement 

U Coefficient cf } 



3.2| Coefficient c t 



(i) 



Non-perturbative tuning ^ 
[4.1| Measurement of c sw 13 
14.21 Non-Perturbative measurement of ca |15 



IHl Conclusions and outlook 17 



1. Introduction 

Quantum field theories with nontrivial infrared fixed points of the f3- function have re- 
cently been studied due to their applications in beyond Standard Model model building. 
In these theories the coupling runs when probed at very short distances, but becomes a 
constant over some energy range in the infrared and the theory appears conformal. One of 
the phenomenological connections is the unparticle [jl], 0, EJ, i.e. the possibility of a fully 
conformal sector coupled only weakly to the Standard Model through effective operators at 
low energies. Another phenomenological motivation to study theories which either feature 
an infrared fixed point or are, in theory space, close to one which does, originates from tech- 
nicolor (TC) and the associated extended technicolor (ETC) models. These models were 
devised to explain the mass patterns of the Standard Model gauge bosons and fundamental 
fermions without the need to introduce a fundamental scalar particle [Q, ||, fj|. 

Early TC models, based on a technicolor sector straightforwardly extrapolated from a 
QCD-like strongly interacting theory, lead to too large flavor changing neutral currents due 
to the extended technicolor interactions. The problems of these simple TC models are solved 
in so called walking technicolor theories || ^, |l0[ |Tl|| . These theories are quasi-conformal, 
i.e. the evolution of the coupling constant is, over a wide range of energy, governed by an 
attractive quasi-stable infrared fixed point at strong coupling. 
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To build walking TC models one needs to tune the gauge and matter degrees of freedom 
so that the desired quasi-conformality arises. To achieve this in SU(iV) gauge theory with 
fermions in the fundamental representation several O(10) Dirac flavors are required. These 
contribute to the precision parameter S, which becomes too large to be compatible with 
the current observations. To obtain enough screening, as required for quasi-conformality, 
but with smaller number of flavors, one considers fermions in higher representations. It has 
been suggested [^] that an ideal candidate for minimal walking technicolor theory would 
be the one with just two (techni)quark flavors in the two-index symmetric representation 
of SU(2) or SU(3). 

Reliable quantitative studies of the models, especially evaluating the ^-functions, re- 
quire lattice Monte Carlo simulations. There are several recent studies of both SU(2) 
13|, [r§ [L5|, [L6|, 0, [E|, fL|, |2|, [21], 



24, 25, 



gauge theories with two- 



and SU(3) || 

index symmetric representation fermions. For related studies in QCD-like theories with 
fundamental representation fermions see p8|, Bi 



, 31 



a m m m 



In this paper we consider the case of SU(2) gauge fields with two fermions in the 
two-index symmetric representation, which, for SU(2), is equivalent to the adjoint repre- 
sentation. So far the lattice studies of this theory have been performed using unimproved 
Wilson fermion action and are hence subject to large 0{a) lattice artifacts. In this paper 
we present the computation of C(a)-improvement. This is a generalisation of the program 
used earlier to compute the improved action for two fundamental representation fermions 
in SU(3) gauge theory [^j], 43, 44, [45|, 46|. The early results of this calculation have 
been presented in refs. |^7|, [Hfl . 

The Wilson fermion action can be improved for on-shell quantities by adding the well- 
known clover term. We tune the coefficient of the clover term (Sheikholeslami-Wohlert 
coefficient ^9|) non-perturbatively, using the Schrodinger functional method. For the 
measurement of the coupling constant we also need the improvement coefficients of certain 
boundary terms. This computation is done using perturbative analysis. For comparison, 
we also calculate the improvement for SU(2) gauge theory with two flavors of fundamental 
representation fermions. 1 

The paper is structured so that in section ^ we first recall the basics of the model as 
well as of the lattice formulation we use. In section |3| we present our perturbative results 
for the boundary terms and nonperturbative results for the improvement coefficients are 
presented in section ^| In section || we conclude and outline the directions of our future 
work. 



2. Lattice formulation: the model and O(a) improvement 

We study SU(2) gauge theory with two different matter contents: two mass-degenerate 
flavors of Dirac fermions either in the adjoint or in the fundamental representation. The 

1 Non-perturbative improvement of the clover term has been recently published for SU(3) gauge field 
theory with 2- index symmetric (sextet) fermions, using the HYP-smeared link clover action ]Hc|]. 
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continuum theory in Euclidean spacetime is defined by 

C=-TcF tu ,F tu , + ^2^ a (Hp + m)il3 a (2.1) 

where F^ is the usual SU(2) field strength, and the gauge covariant derivative is 

D^=(d tl -igA-T a )^ (2.2) 

where a = 1,2,3 and the generators T a are taken either in the fundamental (T a = cr a /2) 
or in the adjoint representation ([T a ] bc = —ie ahc ). The summation in Eq. (|2.1| ) is over 
a = u, d. 

Our main goal in this work is to establish nonperturbative 0{a) improved lattice 
implementation of these theories. While the improvement has been discussed in detail in 
existing literature for SU(3) gauge field with fundamental fermions, the studies of adjoint 
flavors require some alterations. Hence we find it necessary and useful to repeat essential 
parts of the analysis in detail here. 

First recall the usual 0(a) improvement obtained by Sheikholeslami and Wohlert | }49| |. 
The lattice action, split to the gauge and fermionic parts Sq and Sf, is 

So = S G + S F . (2.3) 

Here we use the standard Wilson plaquette gauge action 

Sg = PlY1 (l-^TriO (2.4) 

where /3l = ^/9o an d the plaquette is written in terms of the SU(2) fundamental rep- 
resentation link matrices U^(x), which act as parallel transporters between sites x and 
x + ap: 

Fx\fiv — 

U^x)U v {x + afj,)Ul(x + ai>)Ul(x). (2.5) 

The Wilson fermion action, Sf, for Nf (degenerate) Dirac fermions in the fundamental or 
adjoint representation of the gauge group is 

S F = a A ^2Y,^(x)(iD + m qfi l)Mx), (2.6) 

a x 

where the usual Wilson-Dirac operator is 



D = Mv; + v,)-av;v,), (2.7) 



involving the gauge covariant lattice derivatives V M and V* defined as 

V^(x) = -fyxMx + ap) - ifi(x)], (2. 
a 

V*^(x) = -^{x) - U~\x - ap,)i>(x - ap,)}. (2. 
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Here, the link variables are the usual ones, U^{x) = Un(x), for fermions in the fundamental 
representation while for the adjoint representation they are 

Uf(x) = 2Tr (T a U^x)T b Ul(x)), (2.10) 

where T a , a = 1,2,3, are the generators of the fundamental representation, normalised as 
r £iT a T b = 7}5 ab . We note that in the adjoint representation the elements of [/-matrices are 
real and U^ 1 = U T . 

The lattice action Q2.3| ) is parametrised with two dimensionless parameters, Pl = 
4/g^ are and k = l/[8 + 2am 9) o]. The parameter k is related to the fermion mass. In 
the continuum limit a 4 ~> f d x as a — > 0, and the leading order contribution from 
Q2.3| ) yields the continuum action while the terms of higher order in a will be suppressed; 
these terms are generically termed "lattice artifacts". Since gauge invariance forbids any 
contribution from dimension five operators to the gauge action, only the fermion action 
here is subject to lattice artifacts at 0(a). These are removed (for on-shell quantities) by 
considering the improved action 

5impr = So + 5S SW , (2-11) 
5S SW = a 5 ^2c sw i)(x)^a tlu F tJ , u (x)ij(x) (2.12) 

X 

and tuning the Sheikholeslami-Wohlert coefficient c sw at each (3^ so that the 0(a) effects 
in on-shell quantities cancel; to lowest order in perturbation theory c sw = 1 p9| . Here 
a ihv = ^bVt>7^]/2 and F^ v {x) is the "clover term", lattice field strength tensor in the 
appropriate representation symmetrized over the four /x, ^-plane plaquettes which include 
the point x. 

Because our aim in future work is to measure the evolution of the gauge coupling 
constant using the Schrodinger functional method, we also need to consider the improve- 
ment of the action at the special Schrodinger functional boundary conditions. Schrodinger 
functional method is also used in this work for measuring c S w? but for this the boundary 
improvement is not necessary. 

We consider a system of size L 3 x T, with periodic boundary conditions to the spatial 
directions and with Dirichlet boundary conditions for the gauge fields to the time direction: 

U k (x = 0) = W(k), U k (x = T) = W'(k), (2.13) 

where k = 1,2,3; the explicit form of the boundary fields will be discussed later. For the 
measurement of the coupling constant the boundary gauge fields are chosen so that they 
lead to a constant background chromoelectric field. Due to the frozen boundaries there 
now exists 0(a) contribution to the gauge part of the action, and to account for these we 
consider 

ScM P r = tj- 5>(p)tr(l - U(p)), (2.14) 
v 

where the weights w(p) are equal to 1 for plaquettes in the bulk, w(p) = c s /2 for spatial 
plaquettes at xq = and T and w(p) = ct for time-like plaquettes attached to a boundary 
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plane. The parameters c s and q are tuned to reduce the O(a) boundary contributions. 2 
To leading order in perturbation theory cj = c s = 1. For the electric background field 
which we consider the terms proportional to c s do not contribute. 
The boundary values of the fermion fields are set as 

P+i/>(x = 0,x) = p(x), P-^Oo = r,x) +p'(x), 

P_V>(x = 0, x) = P + i;(x = T, x) = 0, 1 " ' 

with similar definitions on the conjugate fields. The projection operators are P± = |(l±7o). 
The boundary fields p, p' are source fields for correlation functions, and they are set to zero 
when generating configurations in simulations. In the spatial directions it is customary to 
introduce a "twist" for the phase of the fermion fields |S] : 

4){x + Lk) = e idk ip(x), 4>(x + Lk) = 4>(x)e~ tek . (2.16) 

In this work we use 9k = tt/5 throughout. The twist, together with the Dirichlet boundary 
conditions, regulates the fermion matrix so that simulations at zero fermion masses become 
possible. 

The improved lattice action is now given by 

Simpr = SG.impr + Sf + SS SW + SSp,b- (2-17) 

Now the Sheikholeslami-Wohlert term only accounts for the bulk, 

T-a 

5S CW = a 5 ^2 y~] c sw i/j(x)ja l j iV F l j iU (x)^(x), (2.18) 

XQ=a x 

while the boundary effects are captured by SSp,b- This counterterm has two contributions, 
controlled by parameters denoted by c s and q. The term proportional to c s is 



SS, = a 4 (c, - 1) V 



S -j / 



i$(0,x)P_7 fc (V£ + V*)P+V(0,x) 



+^(L,x)P +7fe (V^ + V fe )P_^(L,x) 



(2.19) 



and it clearly vanishes if we set fermionic fields to zero on the boundaries. 

So, similarly to the gauge action, only the term proportional to q contributes, and 
this contribution is given by 

5S Fh = a 4 V(q - l)-^(x)^(x)(5(x - a) + 5(x — (L — a)). (2.20) 
^— ' a 

X 

This can be seen as a correction to the bare mass term at xq = a and xq = L — a, hence 
accounted for by the modification 

mqfl ^ m q ,0 + (q - l){S t , a + $t,L-a)- (2.21) 



2 Recall that gauge invariance guarantees that there are no O(a) contributions to the gauge action in the 
bulk, and hence the boundary terms controlled by c s and Ct are the only ones which arise to O(a) in the 
gauge action. 
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It is known that q = 1 to leading order. 3 

Hence, to obtain 0(a) improvement we need to determine the parameters q, q and 
c sw in the action ( 2.17|) , The parameters Cf and q are determined perturbatively as will 



be described in the following section. The parameter c sw is determined nonperturbatively, 
and this will be determined in section |3|. 

3. Perturbative analysis of the boundary improvement 

As explained in the previous section, due to the Dirichlet boundary conditions associated 
with the Schrodinger functional formalism, we are led to counteract 0{a) lattice artifacts 
on the boundaries both in the gauge and fermion parts of the action. In this section 
we describe in detail the analysis of the rquired counterterms. Although we are mostly 
interested in matter fields in fundamental or adjoint representation of SU(2) gauge group, 
we will present the results applicable also for higher representations of SU(3) since these 
are relevant for the current developments in the studies of these theories on the lattice. 

In principle there exists four counterterms associated with the spatial links in the 
boundary and with temporal links connected to the boundary. Due to the specific form of 
the background field we have chosen, only two of these are needed and these are denoted 
by ct and c±. These boundary coefficients have a perturbative expansion of the form 

c x = \ + c^gl + 0{gi). (3.1) 

Our goal is to determine cj and ct to one-loop order in perturbation theory. 

3.1 Coefficient c| 1} 



We follow the analysis performed in [44] for the fundamental representation. The result of 



|4| is 



c t (1) = -0.0135(1)C F , (3.2) 

and this generalizes to other fermion representations simply by replacing the fundamental 
representation Casimir operator Cf with Casimir operator Cr of the representation R 
under consideration. This is so because the relevant correlations functions are proportional 
to the diagrams presented in figure |l|, which all include the color factor ^2 a {T a ) 2 = Cr. 
Thus it can be shown that also cx Cr. 



In the case of fundamental fermions the original result of [44] is directly applicable with 
Cf = — 1)/(2N C ) = 3/4 for N c = 2. For the other case we have fermions transforming 
in the adjoint representation of SU(2), for which the Casimir invariant is Ca = 2. The 
results for different gauge groups and fermion representations are shown in table [j]. 

3.2 Coefficient cj 1} 

The coefficient cj 1 ^ can be split into gauge and fermionic parts 

&=$»+$*N f . (3.3) 



3 This is so because free Wilson fermions are not subject to 0(a) artifacts. 
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Figure 1: Diagrams contributing to the calculation of The shaded blob on each diagram 

indicates the insertion of the operator T x = {1,75}. 

The contribution c^ 1 ' ^ is entirely due to gauge fields and has been evaluated in [51] for 



SU(2) and in [52] for SU(3). The fermionic contribution cj 1 ' 1 ^ to q has been evaluated 



for fundamental fermions in both for SU(2) and SU(3). We have extended these 
computations for SU(2) and SU(3) gauge theory with higher representation fermions and 
for SU(4) gauge theory with fundamental representation fermions. 

The method we have used is the same as the one presented in |jl3|l , with two exceptions. 
First, the boundary fields have to be transformed to the desired fermion representation. 
Generally the boundary fields are of the form 

U(x, k)\( XQ=0) = exp(aC fc ), U(x, k)\( XQ=L) = exp(aC' k ), (3.4) 

where 

C k = -diag(0i, ... ,4> n ), C k = -diag(0 1; ...,0 n ), (3.5) 

and n is the dimension of the representation. The transformed boundary fields are obtained 
from the fundamental representation counterparts for adjoint representation via ( 2.10| ). 



After the transformation one simply diagonalizes the resulting matrices and ends up with 
a matrix of the form 

diag (exp[i(j)f], . . . , exp[i^]) , (3.6) 
where 4>f give the adjoint representation boundary fields 

C fc A = -^diag«...,^), ^ = -^diag(^,...,0- (3-7) 

For the symmetric representation the components of the boundary fields and C£ 
can be obtained by taking all the symmetric combinations of fa. For SU(3) sextet repre- 
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sentation this is 

4>f = 0i + 0i, 

02 = 01 + 02, 

03=01 + 03, ^ 
04 = 02 + 02, 
0f = 02 + 03, 
06 = 03 + 03- 

The other crucial note comes from the normalization in the calculation of c^ 1 ' 1 ^. Using 
the Schrodinger functional scheme and taking the lattice action with constant background 
field as an effective action r , the running coupling is defined via 

f = f?- 0.9, 

The boundary fields Ck and C' k are functions of the parameter ij so the running coupling is 
given by the change of the system as the boundary fields are altered. The effective action 
T is to one loop order in perturbation theory 

r = 5o" 2 r + ri + o( 5o 2 ), (3.10) 

so the running coupling can be written, as a function of the bare coupling go, i n the form 

=;2 _ „2 / -1 2^;^i 



S=9 2 o{l-9 2 o^r)+0(9 b ). (3.11) 

On small lattice spacings a, the one loop correction Ti diverges. This leads to renormal- 
ization of the lattice coupling, which is given in terms of the bare coupling as 



where z\ = 26q ln(a^) and 



gl t = g 2 + z ig 4 +O(g 6 ), (3.12) 



^C A - \t{R)N f ) (3.13) 



(4vr) 2 V 3 A 3 



is the coefficient in one loop beta function. Now we can write the running coupling as a 
function of the renormalized coupling 



-2 2 
9 =5lat 



a„ri 



+ <Wat)- (3-14) 



The one loop correction to the effective action T± can also be written as 

1 1 

T 1 = - In det Ai - In det A - - In det A 2 , (3.15) 

where the operators Ao and Ai are related to the gauge fixing and pure gauge part of 
the action and the operator A2 = [(D sw + mo)75] 2 is related to the fermionic part of the 
action. The operator D sw is the lattice Dirac operator that includes the Sheikholeslami- 
Wohlert term. Now for the calculation of A2 one needs to transform the boundary fields 
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to the appropriate representation. However in the calculation of Tq one needs to keep the 
boundary fields in the fundamental representation. This is so becouse the pure gauge part 
of c!jp should be independent of the representation of the fermions and this can only be 
achived if the boundary fields in Tq are kept in the fundamental representation. Also this 
produces the expected behavior for the series expansion of 



d v T 



+ zi 



(3.16) 



With these remarks, the numerical calculation is straightforward. The results for 
the nonzero improvement coefficients are tabulated in table |l|. The numbers beyond the 
fundamental representation are new, while those for the fundamental representation provide 
a good check on our computations. For the application to minimal walking technicolor, 
the relevant numbers are the ones on the second line of table |l[ 

Our results are consistent with the generic formula 



0.019141(2T(7?)), 



where T(R) is the normalization of the representation R, defined as Tr(T^T^) 



(3.17) 
T(R)6 ab . 

For the details of the numerical method used to determine coefficient cj: , we refer to the 
original literature where the method was developed and applied first for the pure gauge 
theory case in |5l[] , and later for fundamental representation fermions in [0, [43| . 



N c 


rep. 


Ji.o) 
c t 




~(D 
c t 


2 


2 


-0.0543(5) 


0.0192(2) 


-0.0101(3) 


2 


3 


-0.0543(5) 


0.075(1) 


-0.0270(2) 


3 


3 


-0.08900(5) 


0.0192(4) 


-0.0180(1) 


3 


8 


-0.08900(5) 


0.113(1) 


-0.0405(3) 


3 


6 


-0.08900(5) 


0.0946(9) 


-0.0450(3) 


4 


4 




0.0192(5) 


-0.0253(2) 



Table 1: The nonzero improvement coefficients for Schrodinger functional boundary conditions 
with electric background field for various gauge groups and fermion representations. 

,(1,1) 



We have also plotted our results of cl ±,± ' scaled with \/{2T{R)) against ( [3.17 ) in 



figure |2j Although we did not achieve the accuracy of the original work ]43j , our results are 
fully compatible for fundamental representation fermions. The figure also clearly indicates 
that Cj 1,1 ^ scales with 2T(R). 



4. Non-perturbative tuning 

The continuum physics we are interested in corresponds to massless fermions, so we need to 
simulate at zero physical quark mass. With Wilson fermions the bare quark mass is addi- 
tively renormalized, and the zero of the physical quark mass corresponds to tuning the bare 
quark mass to a critical value, mo = m cr . This tuning is done nonperturbatively and al- 
lows for determination of the improvement coefficient c sw simultaneously. Here we describe 
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Figure 2: Our results of c\ scaled with 2T(R) compared with conjectured value of 
4 M) /(2T(i?)). 



the calculation of c sw for Nf = 2 flavors of SU(2) fundamental and adjoint representation 
fermions. 

In these simulations the fermion fields have the boundary conditions given in Eqs. Q2.15IJ2.16D 
For the fundamental representation fermons we fix the gauge field Dirichlet boundary con- 



ditions at xo = and xq = T [41]: 

U k (x = T) = exp(iC"), C 
U k (x = 0) = exp(iC), C = 



-/r aa 
37r acx 3 

'T~zT : 



(4.1) 
(4.2) 



for k = 1,2,3. Because the boundary link matrices commute, we call these boundary 
conditions Abelian, in contrast to the non-Abelian (non-commuting) ones defined below. 

The physical quark mass is defined via the partial conservation of the axial current 
(PCAC) relation, 



where 



, 1 |(gg + do)f A (x ) + c A ad* Q dof P (x ) 

M(x ) = -r-, — t = r(x ) + c A s{x ), 

2 Jp{xo) 



K = ^{x)l^^-a a ^{x), 



P a = ii(x) l5 -a a il>(x) 



f A (x ) = -a 6 ^(^(x)C(y) 7 5^ a C( Z )>, 
f P (x ) = -a 6 ^(P a (x)C(y)75^\(z)). 



(4.3) 

(4.4) 
(4.5) 
(4.6) 

(4.7) 



- 10 - 



Another set of correlation functions, f' A and f' P is denned via 

f' A (T - xo) = -a 6 ^(AS(x)C(y)75^ a C(z)}, (4.8) 

f P (T - x ) = -a 6 ^(P a (x)C(y) 75 ^ a C(z)). (4.9) 

The bare mass is tuned so that M(T/2) vanishes. The c sw term is tuned simultaneously 
using mass measurements at a different point in the bulk looking for variations of the order 
of the lattice spacing. Defining M' with obvious replacements of primes, it follows that the 
quantity 

AM(x ) = M(x ) - M'(x ) (4.10) 

vanishes up to corrections of 0(a 2 ) if both c sw and ca have their proper values. In order 
to recover the correct tree level behaviour we fix these quantities M and AM two their 
tree level values, measured by from a cold gauge configuration with k c = 0.125. This gives 
a small correction to the relations: 

AM(x ) = M(x ) - M'(x ) -5 = 0, M(x ) = 5 M (4.11) 

However, for the adjoint representation fermions there are complications which sig- 
nificantly reduce the effectiveness of the above method. Using Eq. ( [2.10D we immediately 
notice that the Abelian boundary matrices ( |4.1|]4.2j ) are transformed into form 

/ o\ 

U k = (4.12) 

V V 

Thus, there is a component of the adjoint representation color vector which simply does 
not see the background field. This feature is independent of the color structure chosen 
for the boundary conditions. It turns out that regardless of how the fermion sources or 
the constant boundary conditions are chosen, at long distances the correlation functions 
behave as if there is no background field. In other words, the adjoint fermion correlation 
functions "see" the background electric field only at short distances. This significantly 
reduces the effectiviness of the background field method for tuning c sw . 

This effect can be improved by using boundary conditions which maximize the dif- 
ference between the two boundaries. We use the following asymmetric "non- Abelian" 
boundary conditions: links at the upper xo = T boundary are chosen to be trivial 

U{x = T,k) = I (4.13) 

and at the lower boundary xq = we use 

k 

U(x = 0,k) = eMaC k ), C k = ^- T —. (4.14) 

2 iL 

This creates a strong chromomagnetic field at xq = boundary. These boundary conditions 
do not fully cure the problem, but nevertheless provide enough leverage so that the PCAC 
mass relation can be used to tune c sw . 
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Figure 3: Left: Fundamental representation fermion mass aM(xo) measured from the classical 
gauge field configuration satisfying the Abelian boundary conditions ( |4.1| , 4.2) on a 8 3 x 1 6- lattice. 
Bare mass is amo = 0.01, which is also aM in the continuum limit. Inclusion of the clover term 
(c sw = 1) significantly reduces the cutoff effects. Right: aAI(xo) for adjoint representation fermions 
and for the Abelian boundary conditions (hi), (h2), and for the non- Abelian boundary conditions 
( 4.13 , 4.14 ) . Here the correlation functions between c sw = and c sw = I differ significantly at long 
distances only for the non-Abelian boundary conditions. 



This behaviour can be demonstrated already at the classical level: in figure ^ we show 
the PCAC fermion mass (|4.3|), measured using the classical minimum action gauge field 
configuration which satisfies the appropriate boundary conditions. The bare fermion mass 
has been set to amo = 0.01, and, in the absence of the background field or lattice cutoff 
effects, the PCAC measurement would yield precisely this value. However, with finite 
lattice spacing the non-trivial classical background field gives rise to cutoff effects, which 
moves the PCAC mass away from aM = 0.01. For the fundamental representation fermions 
and the Abelian boundary conditions fl4.1| ), ( |4,2|) . (left panel in figure |3|), we can observe 
that setting c sw = (non-improved standard Wilson fermions) the measured mass values 
are far from the continuum limit, whereas using c sw = 1 (the correct value at the classical 
level) these effects are strongly reduced. 

For the adjoint representation fermions the behaviour is very different, as shown on the 
right panel of fig. Hf using the Abelian boundary conditions the measured masses aM(xo) 
rapidly approach 0.01 /a as xq increases, for both c sw = or 1. This indicates that the 
correlation function lacks the sensitivity to c sw and cannot be used for tuning it to the 
correct value. 

On the other hand, with the non-Abelian boundary conditions (4.13, 4~T4| ) the corre- 
lation function remains sensitive to the value of c sw to longer distances. The sensitivity 
remains in the mass asymmetry AM(xo), ( |4.10 ), which can now be used to tune c sw . We 
note that these boundary conditions are useful only for determining c sw , not for evaluating 
the coupling constant. 
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In order to remove the dependence on ca, for fundamental fermions, we consider 

M(x , y ) = r(x ) - S (x ) T< f°\ - ^ , (4.15) 

which coincides with M(xq) up to 0(a 2 ) corrections and is independent of ca- With adjoint 
fermions this quantity suffers from large statistical fluctuations and is not useful. Instead 
we simply consider the quantity M(xq) and fix ca to its perturbative value [44] 



c A = -0.00567(1)CW + 0(g 4 ). (4.16) 

We then measure ca separately to confirm the validity of our choice. 

In order to evaluate c sw we used the following routine: we choose lattice volume L 3 x 
T = 8 3 x 16 for both fundamental and adjoint representation fermions, and a set of values 
of the lattice coupling (3. For fundamental fermions we measure M = M(T/2,T/4) and 
AM = AM(3T/4,T/4). For adjoint fermions we measure M = M(T/2) and AM = 
AM(3T/4) fixing ca to its perturbative value. 

1. For a given j3, we choose initial c sw (typically extrapolating from results obtained 
with previous values of /?). 

2. We choose a couple of values for k, and determine by interpolation the critical value 
k c (/3, c sw ) where the fermion mass M is equal to the tree level value. 

3. Once we have an estimate of the critical k, we choose a new value for c sw and repeat 
the search of k c . 

4. At the same time, we measure AM(c sw ). Now we can linearly interpolate/extrapolate 
in c sw so that AM vanishes, obtaining the desired value of c sw . Using simulations at 
this final c sw we can relocate the critical k, if desired, and verify the results of the 
interpolation. 

The above tuning is done at small L/a, and the results are applied for all lattice sizes 
since the L/a dependence is expected to be weak. Furthermore, we only consider a range 
of P and fit the critical values to an interpolating function to obtain m c (/3i) and Cg W (/3x). 

4.1 Measurement of c sw 

In figures |I| and [5] we show our results for the clover coefficient c sw for both fundamental 
and adjoint representations. The values of (3 used are f3 = 2.5, 3, 4, 5, 6, 8, and also (5 = 
2.25 and 10 for the adjoint representation. To clarify the tuning method we provide the 
measurements of M and AM with adjoint fermions in table 0. In tables || and |I| we give 
our results for c sw for fundamental and adjoint fermions respectively. 

Finally, the measured values for c sw can be fitted with a rational interpolating expres- 
sion, which can used in simulations for this range of /3-values. For fundamental represen- 
tation fermions we use the perturbative 1-loop result c sw = 1 + 0.1551(l)g 2 + 0(g A ) to 
constrain the fit: 

1 - 0.090254 5 2 - 0.0388465 4 + 0.028054# 6 



1 - (0.1551 + 0.090254)c/ 2 



(4.17) 
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P 

Figure 4: c sw for two flavors of fundamental representation fermions. The solid line is the inter- 



polating fit, Eq. (4.17), and the dashed line is the 1-loop perturbative value 




Figure 5: c sw for two flavors of adjoint representation fermions, with the interpolating fit, 



Eq. (4.1 



For the adjoint representation the perturbative result is not known, and we obtain the fit 
result 

1 + 0.032653# 2 - 0.002844c/ 4 



1 -0.3141535 2 



(4.18) 



In both cases the interpolating fits are valid for j3 J> 2.5. For the adjoint fermions it is 
difficult to reach smaller /3-values because c sw grows rapidly, and while we were able to 
reach /3 = 2.25 the errors were too large to constrain the fit ( f4.18| ) further. 
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Table 2: Results for the quark mass M and AM with two fermions in the adjoint representation 
4.2 Non-Perturbative measurement of ca 

When measuring c sw for adjoint fermions we chose to keep the coefficient ca at the pertur- 
bative value. In figure ^| we show how the choice of ca affect AM(3T/4) with certain choice 
of parameters. Typically ca is between —0.005 and —0.01 at the range of B we explored. 
We see that even differences of this order have small effect to AM. 
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Table 3: Results for c sw with two flavors of fermions in the fundamental representation 
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Table 4: Results for c sw with two flavors of fermions in the adjoint representation 
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Figure 6: The dependence of AM(3T/4) of ca- The measurement was done with /3 = 4 and 
Csw = 1.522. 



To verify the accuracy of our choice we have also estimated a non-perturbative value 
for ca- For this we have used the same improvement condition as in BE]. We do simulations 



with two different values of the fermion phase 6 in the boundary conditions [2.16| , using the 
measured values of c sw and k c above. Without any discretisation errors the difference in 
the measured masses should be equal to the tree level value. Requiring that this condition 
is met, we can find an estimate of ca- 

From two simulations with 8 = and 9 = ir/2 we calculated the discretisation effect 



AM(ca)' = M(x 



0, c A ) -M(x = 8;9 = it/2, c a ) - S, 



(4.19) 



where 5 is the tree level value of the difference. It is similar to the tree level correction in 
equation 4.11 and is relatively small. These simulations were done using a trivial boundary 
condition, where all the boundary matrices were set to unity. Depending on the lattice 
coupling between 2000 and 35000 trajectories were performed for each value of 0. 
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Figure 7: ca for two flavors of adjoint fermions. The solid line is the interpolating fit, Eq. (fo?]), 
and the dashed line is the 1-loop perturbative value 

As the quark mass, and therefore AM' is simply linearly dependent on ca, we can 
measure AM' for two of values of ca to find the correct value where AM' = 0. 
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Table 5: Results for ca 



The results for ca are given in table |5] and depicted in figure [j]. We see that in the 
region where we have measured c sw and ca, it is justified to use the perturbative value for 
ca- 

5. Conclusions and outlook 

We have calculated O-improvement of SU(2) gauge theory with two Wilson fermions in the 
fundamental or adjoint representation. The main results are the non-perturbative evalua- 
tion of the Sheikholeslami-Wohlert clover coefficient c sw and the perturbative calculation 
of the boundary improvement terms needed for full improvement in the Schrodinger func- 
tional formalism. The result for c sw is generally applicable to lattice simulations of these 
theories. We also verified that the axial current improvement coefficient ca is well described 
by the 1-loop perturbative formula in the range of lattice spacings studied. In addition to 
the perturbative results on SU(2) gauge theory and adjoint fermions, we obtained results 
also for SU(3) and adjoint or sextet fermions which will be useful also for other groups 
studying these theories. 
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The main application for the improved action is more accurate lattice Monte Carlo 
analyses of the candidate theory for minimal walking technicolor, SU(2) gauge theory with 
two adjoint representation fermions. The boundary improvement terms permit improved 
measurement of the evolution of the coupling constant with the Schrodinger functional 
scheme. Indeed, in earlier unimproved analyses 17, [l^] significant cutoff effects were ob- 
served at coarse lattices. The measurement of the coupling with the improved action is left 
for future work. 
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